Superfluidity and Quantum Melting of p-li2 Clusters 
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Structural and superfluid properties of P-H2 clusters of size up to A'^=40 molecules, are studied 
at low temperature (0.5 K < T < 4 K) by path integral Monte Carlo simulations. The superfluid 
fraction ps{T) displays an interesting, non-monotonic behavior for 22 < < 30. We interpret this 
dependence in terms of variations with A'^ of the cluster structure. Superfluidity is observed at low T 
in clusters of as many as 27 molecules; in the temperature range considered here, quantum melting 
is observed in some clusters, which are seen to freeze at high temperature. 

PACS numbers: 67.90.+Z, 61.25.Em 



Recent developments in spectroscopy afford the inves- 
tigation of properties of a single complex molecule em- 
bedded in clusters of ''He or P-H2 . Specifically, by study- 
ing the rotational spectrum of the molecule one can ob- 
tain evidence of decoupling of its rotation from the sur- 
rounding medium (i.e., the cluster), at sufficiently low 
temperature (of the order of a fraction of 1 K). Such a 
decoupling is interpreted as due to onset of superfluidity 
(SF) in the cluster 0. This is now an area of intense, 
current research effort, aimed at gaining theoretical un- 
derstanding of the microscopic origin of SF, perhaps the 
most fascinating manifestation of quantum behavior on 
a macroscopic scale. In particular, theoretical questions 
are being addressed such as: What is the smallest fi- 
nite size system for which SF can be observed ? Which 
condensed (i.e. non-gaseous) matter systems, besides he- 
lium, can display this phenomenon, if not in the bulk at 
least in sufficiently small clusters ? 

With respect to the second question, droplets of hydro- 
gen molecules are clearly of fundamental interest. Molec- 
ular para-Hydrogen has long been speculated to be a po- 
tential superfluid, owing to the bosonic character and the 
light mass of its constituents 0]. However, the search for 
SF in bulk P-H2 has so far been frustrated by the fact 
that, unlike helium, this system solidifies at low tem- 
perature, as the intermolecular potential is significantly 
more attractive than that between two helium atoms. On 
the other hand, clusters of P-H2 molecules of sufficiently 
small size, ought to remain "liquidlike" at significantly 
lower temperature than the bulk, possibly turning su- 
perfluid [3]. 

Indeed, theoretical studies carried out some fifteen 
years ago, based on path integral Monte Carlo (PIMC) 
simulations 0, yielded evidence of a finite superfluid re- 
sponse in pure [(p-H2)Ar] clusters, with A^=13 and iV=18 
molecules, at a temperature T < 2 K, whereas a larger 
cluster (A^=33) was found to be "sohdlike", nonsuper- 
fluid, in the same temperature range. Other PIMC calcu- 
lations showed that, while SF may not necessarily occur 
in larger clusters, there is nonetheless a large propensity 
for quantum exchanges Finally, more recent PIMC 
studies have yielded evidence of superfluid behavior in 
small (17 molecules) clusters of P-H2 doped with a single 
OCS 1113 or CO Q molecule. 



For pristine P-H2 clusters, the only theoretical study of 
superfluidity is restricted to three cluster sizes (A^=13, 
18, and 33) 0. No systematic study of the superfluid 
properties of clusters as a function of size and temper- 
ature has yet been carried out. Although experimental 
data are not yet available for pure P-H2 clusters, novel 
techniques based on Raman spectroscopy hold promise 
for the investigation of superfluidity in these systems ^ . 

In this Letter, we present a detailed investigation of (p- 
H2)Ar clusters at low temperature (down to T—Q.5 K), for 
N < 40, by means of PIMC simulations, based on a re- 
cently developed worm algorithm [lol |. This numerical 
technique affords accurate estimates of thermodynamic 
properties of Bose systems. In particular, the super- 
fluid fraction can be calculated with much greater ac- 
curacy than that afforded by conventional PIMC, which 
has been one of the leading many-body computational 
methods of the last 20 years ^ . 

Our main findings are the following: (p-H2)7v clusters 
with N < 22 are liquidlike, and superfiuid at low T. Su- 
perfluid properties of clusters with 22 < N < 30, strongly 
depend on N. A few clusters, in the range 22 < N < 30, 
feature, in temperature range considered here, a behavior 
that is consistent with coexistence of insulating (solidlike) 
and superfluid (liquidlike) "phases" , the latter becoming 
prominent as T is lowered. In other words, such clusters 
melt at low T as a result of zero-point motion, and freeze 
at higher temperature. We refer to this intriguing be- 
havior as "quantum melting" . The superfluid response 
of clusters with > 30 is significantly depressed; no- 
tably, however, for = 40 permutation cycles can still 
be observed including a considerable number (as many 
as 20) of P-H2 molecules. 

We model our system of interest as a collection of A^ 
P-H2 molecules, regarded as point particles and inter- 
acting via an accepted pair potential The system 
is assumed to be at a temperature T=l/(3. Since we 
are interested in studying properties of the clusters as a 
function of A^, we use a variant of the worm algorithm 
described in Ref. [lol |. in which the number of particles 
in the configurations inside the so-called Z sector (those 
that contribute to the expectation values of physical ob- 
servables) is fixed at A^. We utilized a high-temperature 
approximation for the many-body density matrix accu- 
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FIG. 1: (color online) Superfluid fraction ps{T) for clusters 
of 20 (circles) and 23 (boxes) P-H2 molecules. Dotted lines 
are guides to the eye. When not shown, statistical errors are 
smaller than the symbol size. 
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FIG. 2: (color online) Superfluid fraction versus cluster size 
A^, at T—1 K (flUed circles). When not shown, statistical er- 
rors are of the order of, or smaller than the symbol size. Solid 
line is only a guide to the eye. Also shown for comparison are 
results from Ref. (open triangles). 



rate up to fourth order in the imaginary time step e [13 , 
and used for our calculations a value of the imaginary 
time step e = 1/640 K^^, which we empirically found 
to yield converged estimates. We computed cluster ener- 
getics, radial density profiles and the superfluid fraction 
Ps{T) (using the well-known "area" estimator ^j). 

Figure n shows ps{T) for clusters of size iV=20 and 
A^=23. In both cases, as expected ps{T) is a monotoni- 
cally decreasing function of T. As we discuss below and 
later in the manuscript, however, the physical behavior 
of these two clusters is qualitatively different. 

For iV=20, the behavior of ps(T') is close to that ob- 
served in Ref. y for N=18. At T < 1.25 K, the system is 
essentially entirely superfluid; ps (T) drops rather quickly 
to a value ~ 0.2 at T=2.5 K (corresponding to roughly 
four molecules in the superffuid phase), and decreases 
more slowly at higher temperatures. While there is ob- 
viously no real phase transition in a finite system, the 
notion of "superfluid fraction" becomes scarcely mean- 
ingful, when the average number of molecules in the su- 
perfluid phase is of the order of 1. Therefore, we opera- 
tivcly define our "transition temperature" Tc as that at 
which Nps{Tc) ~ 2. For iV=20, this heuristic criterion 
yields Tc ~ 3 K. Remarkably, however, even at this tem- 
perature the probability for a P-H2 molecule to belong 
to a permutation cycle involving three or more molecules 
is still as large as ~ 3%, and exchange cycles involving 
as many as 13 molecules are observed. This is consistent 
with the qualitative observation made in Ref. '3, based 
on a PIMC simulation which did not explicitly include 
exchanges. 

As the number N of particles increases, the physics of 
a cluster ought to approach that of the bulk; whereas this 
means a liquid (superfluid at low T) for helium clusters, 
P-H2 forms an insulating crystal. It is therefore reason- 
able to expect that, at some fixed, low T (e.g., T=l K), 
the superfluid fraction of P-II2 clusters should decay to 
zero as ^ 00. 



Figure |21 shows our estimates for the superfluid frac- 
tion PS as a function of the cluster size A'^, at T—1 K. For 
A^ < 22, (p-H2)a' clusters are entirely superfluid (within 
the precision of our calculation) i.e., ps ~ 1- Further- 
more, the temperature dependence of ps is consistently 
similar to that shown in Fig. ^(for A^=20). In the range 
22 < A^ < 30 a remarkable, nonmonotonic behavior of ps 
is observed, with dramatic differences between clusters 
differing by the addition of just one molecule. We pro- 
pose that such an interesting behavior reflects alternating 
liquidlike (superfluid) or solidlike (insulating) character 
of the clusters, strongly dependent on N. This interpre- 
tation is consistent with the expectation that a crystalline 
(nonsuperfluid) phase should emerge, at large A^. 

For example, at A^=25 the cluster displays a large su- 
perfluid response at T~l K. If a single molecule is added, 
the superfluid fraction drops abruptly, to a value less than 

0. 1. On adding yet another molecule, the superfluid frac- 
tion grows again, to approximately 0.25. It seems difficult 
to imagine that the addition of a single molecule would 
alter so drastically the superfluid component, if the clus- 
ter structure stayed essentially the same, e.g., liquidlike. 
It seems reasonable, instead, to relate the changes in the 
superfluid properties to structural changes that occur on 
adding molecules. 

In order to illustrate this point, we show in Fig. 13 
proflles of radial density {p{r)), computed with respect 
to the center of mass of the cluster, at T=l K. For A^=15, 
the large value of p(r — > 0) and the local minimum for r « 
2 A indicate the presence of a single P-H2 molecule in the 
center of the cluster. Other molecules form a floppy shell 
around the central one, as shown by the peak at r—A A. 
On increasing the cluster size, qualitative changes occur 
at A^ ~ 22. The value of p{r — > 0) becomes negligible, 

1. e., the center of the cluster is no longer occupied by a 
molecule. There is a peak at about 2 A from the center, 
as an inner molecular shell forms. A second, broader peak 
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FIG. 3; (color online) Radial density computed with respect 
to the center of mass, for clusters with 15, 25, 26, and 27 P-H2 
molecules. Statistical errors, not shown for clarity, are of the 
order of 5xlO~^ A~"^ or less. 

at larger distance (r « 5 A) corresponds to the formation 
of an outer shell. 

The main structural change, going from N=2b to 
is that the first peak becomes significantly 
sharper, and its height increases by some 40% (see Fig. 
O. We interpret this as evidence that the inner shell be- 
comes more solidlike, with molecules localized and quan- 
tum exchanges depressed, both in the first shell as well as 
between the first and second shells. If another molecule is 
added, the density profile for N=27 features a first-shell 
peak and an intershcU minimum of heights intermediate 
between those of the N=25 and N=26 cases, and ps in- 
creases to a value much lower than for N=25, but signif- 
icantly greater than that for N=26. Thus, the addition 
of a molecule to the N=26 has the effect of frustrating 
the solid order of the inner shell, increasing molecule de- 
localization and leading to quantum exchanges. 

A particularly intriguing behavior is observed in some 
clusters; we discuss in detail the A^=23 one, for which 
Ps{N) at r=l K takes on a local minimum. Fig. 0| 
shows the values of the superfluid fraction (upper panel) 
as well as of the potential energy per particle {V, lower 
panel), recorded in a typical Monte Carlo run (results 
shown in the figure correspond to a few hundred hours 
of CPU time on a high-end workstation). Specifically, 
data shown refer to successive block averages of ps and 
V, each block consisting of 500 sweeps through the sys- 
tem The quantity ps displays large fluctuations, but 
generally switches abruptly between two clearly identifi- 
able regimes, namely, one in which it oscillates around an 
average value of one, and another in which it stays close 
to zero. In close correspondence with ps, the potential 
energy V displays similar oscillations, switching between 
configurations of lower (for small ps) and higher (greater 
Ps) values. 

We interpret this pattern as the signature of the coex- 
istence of two cluster phases, characterized by large and 
small superfluid response. The behavior of the potential 
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FIG. 4: (color online) Behavior of superfluid fraction (up- 
per panel) and potential energy per molecule (lower panel) 
observed during a typical Monte Carlo run for a cluster of 
N—23 molecules at T—1 K. Data shown refer to successive 
block averages of ps and V. A single "block" consists of 500 
sweeps through the entire system (see text). Although large 
fluctuations are present, visual identification is relatively easy 
of two different regimes, in which ps is on average either close 
to 1 or zero; correspondingly V oscillates between two values 
close to ~ —55 K (liquid) and ~ —60 K (solid) respectively. 




FIG. 5: (color online) Radial density profiles computed for a 
cluster of N=23 molecules at T=2.0 K and T=0.75 K. Statis- 
tical errors, not shown for clarity, are of the order of 5xl0~* 
or less. 

energy suggests that, in its non-superfluid phase, a clus- 
ter comprising 23 P-H2 molecules should feature distinct 
solidlike properties, chiefly a high degree of localization of 
the molecules. The coexistence of these two (solid- and 
liquidlike) phases renders the precise determination of 
the average value of ps computationally rather demand- 
ing, i.e., fairly lengthy runs are needed. On decreasing 
the temperature, the liquidlike superfluid phase becomes 
dominant, i.e., the cluster "melts" at low T due to quan- 
tum zero-point motion of the P-H2 molecules, "freezing" 
instead at higher temperature. This is consistent with 
the observed evolution of the radial density profile, shown 
in Fig. for T=2.0 K and 0.75 K. As T is lowered. 
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the first peak broadens significantly, as molecules enjoy 
greater mobility. 

Qualitatively similar results are seen for other clus- 
ters, e.g., N—27, in the temperature range explored in 
this work. It is likely that other clusters, in the range 
22 < < 30, may display the same behavior, at some 
temperature (in some cases possibly much lower than 
r=0.5 K, which is the lowest considered here). It is im- 
portant to note that this behavior is markedly different 
than that observed in clusters with N < 22, for which a 
plot such as that of Fig. 0] merely shows the two quan- 
tities PS and V fluctuate around their average values, 
with no evidence of the system switching back and forth 
between two distinct phases. Indeed, clusters with N< 
22 are found to be liquidlike at all temperatures, with a 



growing superfluid (normal) component at low (high) T. 

Summarizing, we have studied superfluid and struc- 
tural properties of P-H2 clusters of size < 40. We 
observed nontrivial superfluid behavior of (p-Il2)jv, as a 
function of N. Our observation is consistent with the 
emergence of a solid phase, as the size of the cluster 
grows; however, this occurs non-monotonically. Some 
clusters (e.g., A^=23 and 27) feature, at low tempera- 
ture, quantum melting, induced by zero-point motion; 
these clusters are observed to freeze at high temperature. 
Some of these predictions may soon be tested experimen- 
tally gj. 
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